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Preface 


The  woik  described  herein  was  conducted  by  Dr.  Raymond  Chapman  of 
Ray  Chapman  and  Associates.  Vicksburg.  MS.  for  the  U.S.  Army  ^gineer 
Waterways  Experiment  Station  (WES),  under  Contract  No.  DACW39-92-C-(X)- 
14.  Mr.  Teny  K.  Gerald  of  AScl.  Inc.,  developed  the  interactive  interface 
software.  The  contract  was  monitored  by  Dr.  Marie  S.  I>ortch.  (Thief.  Water 
(Quality  and  Contaminant  Modeling  Branch.  Environmental  Laboratory  (EL). 
WES.  Suiq)Ott  for  this  wotk  was  provided  by  the  U.S.  Army  Engineer  Dis¬ 
trict,  New  York,  under  the  Water  Resource  I^velopment  Act  of  1986  (Public 
Law  96-662.  Section  728a). 

This  report  was  prepared  by  I^.  Chapman.  Mr.  Gerald,  and  Dr.  Dortch 
under  the  general  supervision  of  Mr.  Dc^d  L.  Robey.  (Thief.  Environmental 
Processes  and  Effects  Division.  EL.  and  Dr.  John  Harrison,  Director,  EL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert 
W.  Whalia  Commander  was  COL  Bmce  K.  Howard,  EN. 


1  introduction 


Within  the  scope  of  .the  Water  Resource  Development  Act  of  1986  (PL  96- 
662,  Sec.  728a),  a  feasibility  smdy  was  authorized  for  the  purpose  of  develop¬ 
ing  a  monitoring  and  modeling  strategy  for  water  quality  management  of  tte 
New  York  Bight  As  pan  of  the  overall  modeling  strategy,  a  numerical  model 
c^>able  of  predicting  the  transpon  and  fate  of  floatables  and  suspended  or 
dissolved  materials  using  a  thm-dimensional  (3D)  particle  tracing  method¬ 
ology  was  developed.  Specific  features  implonent^  include  (a)  surface  trans¬ 
pon  of  floatables  by  the  wind,  and  (b)  3D  advection  and  diffusion  of 
suspended  and  dissolved  matter.  The  particle  tracking  model  is  interfaced  with 
the  U.S.  Army  Engineer  Waterways  &(periment  Station’s  3D  hydrodynamic 
model  (CH3D)  (Johnson  et  al.  1991).  In  addition,  interactive  window-based 
software  has  been  developed  to  assist  in  the  setup  and  operation  of  the  stand¬ 
alone  version  of  the  particle  tiaddng  model  This  repon  presents  the  particle 
tiatking  model  develo{Hnent  and  verification.  It  is  one  of  five  reports  of  the 
New  York  Bight  Study. 
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2  Model  Formulation 


Panicle  tracking  within  the  CH3D  boundary-fitted  mode^  is  performed  on 
the  transformed  grid  to  take  advantage  of  the  unit  square  computational  cells. 
Particle  positions  x,  are  updated  from  time  level  to  ^  -f  /  as  follows; 


Xi 


Xi  *  Aii 


N 


(1) 


where  the  transformed  contravariant  particle  displacement  vectors  are  approxi 
mated  via  the  following  second-order  Taylor  series: 


Ax. 


«,Ar  + 


dUiAf 
dt  2 


(2) 


where 

Ax,  =  displacement  vector 
u,  =  contravariant  velocity  vet^r 


r  =  time 


and 


dU;  du.  du. 

W  ^  ~3r  * 


(3) 


At  the  end  of  a  desired  ouqxit  interval,  the  contravariant  position  vectors  ate 
transformed  to  Canesian  coordinates  and  stored  altmg  with  grid  ceU  indices  for 
graphical  display  purposes. 

Specific  issues  addressed  during  the  development  of  a  CH3D-based  panicle 
trad^g  model  were  boundary  ctmditions,  the  required  accuracy  of  die  spatial 
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interpolation  and  temporal  integration  schemes,  wind-generated  surface  current 
and  wave  drift,  model  verification,  and  the  mode  of  operation  of  the  particle 
tracking  routine. 

The  required  accuracy  of  the  spadal  interpolation  and  temporal  integration 
schemes  were  examined  by  testing  alternative  versions  of  the  particle  tracking 
model.  The  interpolation  of  the  velocity  components  within  the  3D  grid  is 
accomplished  by  either  a  bilinear  or  biquadratic  interpolation  on  each  horizon¬ 
tal  grid  c&U  and  a  two  point  interpolation  between  vertical  layers  (Hildebrand 
19S6).  Shift  operators  are  employed  to  properly  orient  the  unit  square  interpo¬ 
lation  functions  on  the  space-staggered  hydrodynamic  grid.  The  temporal 
integration  options  for  Equation  2  are  a  ftrst-order  Euler  and  a  second-order 
trapezoidal  (Henrici  1964).  When  the  second-order  time  integration  options 
are  employed,  an  additional  option  is  provided  to  either  include  or  omit  the 
CH3D  convective  acceleration  veaor  (F1,FJ)  in  computation  of  the  particle 
trajectory  (Equation  3). 

Depending  on  the  selection  of  the  vertical  grid  spacing,  the  near-surface 
velocities  predicted  by  CH3D  can  be  several  meters  below  the  free  surface. 

As  a  consequence,  an  Ekman  surface  current  extrapolation  (Neumarm  and 
Pierson  1966)  is  required  to  estimate  the  surftux  drift  requited  for  floatables. 
The  extrrqxrlation  technique  is  based  on  Ekman  surface  drift  velocity  profiles 
which  read  as  follows: 


KU)  *  Vje^COS(i^2)  (4) 

and 

V,(z)  =  V,e^S/iV(pz)  (5) 


where 

V,  =  X  velocity  component 
V,  =  surface  drift  velocity 
a  =  decay  factor  widi  depth  z 
P  =  rotatian  factor  with  depth  z 
Vy=Y  velocity  component 
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Specifically,  the  magnitude  and  direction  of  the  surface  velocity  is 
computed  by  sim{dy  fitting  the  decay  and  rotation  factors  using  the  velocity 
components  (V,  .Vy  defined  within  the  top  two  CH3D  grid  cells. 

The  wind-forcing  methodology  adopted  for  floaud)les  is  based  on  a  balance 
of  die  aerodynamic  and  hydrodynamic  forces  exerted  on  a  floating  object  and 
the  object’s  acceleration.  The  balance  equation  for  the  X  component  of  parti¬ 
cle  motion  is  as  follows: 


at/, 

-IT 


=  /V  *  F  +  F 

^  p  m  w 


(6) 


where 

Ui,  =  X  component  of  the  wind-driven  particle  velocity 
/=  Coriolis  parameter 

V,  =  y  component  of  die  wind-driven  particle  velocity 

F,  =  combined  shear  and  form  drag  force  of  the  wind 

F,  =  comlxned  shear  and  form  drag  force  of  the  water 

A  quadratic  drag  law  is  adopted  for  the  contribution  of  shear  and  form  drag 
forces.  As  a  result,  wind  forces  per  unit  mass  are  written  as  follows: 


where, 

=  air  and  paitide  denaty 
=  skin  fdcdon  and  form  drag  coefficients 
H  =  nominal  particle  height 
L  =  nominal  particle  lengdi 
=  relative  wind  speed,  u.-u, 

Q,  =  resultant  relative  wind  speed 
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An  identical  force  fonnuladon  is  used  for  the  water  forces  with  the  excep- 
ti<Ht  that  water  and  particle  densities  are  used  and  the  sign  of  the  relative 
particle  speed  is  changed.  Using  the  surface  velocities  estimated  by  Ekman 
extrapolation  subroutine,  Etpiation  6  is  integrated  to  obtain  velocity  used  to 
update  the  X  trajectory  disiAacement  equation  (first  term  tm  right  side  of  Equa¬ 
tion  2)  in  the  particle  tracking  subroutine.  The  second  term  on  right  side  of 
Equation  2  is  not  used  for  floatables.  A  similar  equation  for  the  Y  component 
velocity  has  also  been  implemented.  In  addition  to  the  direct  wind  force  bal¬ 
ance,  the  influence  of  wave-induced  transport  (i.e.,  stokes  drift)  is  included  by 
adding  1.5  percent  of  the  wind  speed  to  the  surface  particle  velocity  (Kenyon 
1969). 

Subroutine  WFORCE,  which  estimates  the  wind-driven  transport,  is  called 
by  the  particle  tracking  subroutine  each  time  step.  The  switches  IFLT,  ISURF, 
and  IWIND  are  set  to  one  to  enable  the  call  to  WFORCE.  Data  that  are  re»i 
into  the  particle  tracking  routiiM  are  the  skin  fhctiai  and  form  drag  coeffi¬ 
cients  for  air  and  water,  the  nominal  height  atKl  length  of  the  particle,  and  the 
fraction  of  the  particle  height  that  is  exposed  to  direct  wind  forcing. 

With  regard  to  particle  settling,  a  single  settling  velocity  representative  of 
the  mean  particle  size  or  density  is  specified.  Subsequent  to  nondimensional 
scaling,  the  constant  settling  sp^  is  simply  subtract^  from  the  interpolated 
vertical  velocity  prior  to  the  particle  position  update.  At  tlte  present  time, 
when  a  particle  reaches  the  bottom,  it  remains  there. 

Particle  diffusion  is  assumed  to  follow  a  simple  randmn  walk  process 
(Ftsdter  et  al.  1979).  Specifically,  a  diffusion  distance  defined  as  the  square 
root  of  the  product  of  an  input  diffusion  coefficient  and  the  time  step  is 
decomposed  into  X  and  Y  displacements  via  a  random  direction  function.  The 
Z  diffusion  distance  is  scaled  by  a  random  positive  or  negative  direction.  The 
equations  for  the  horizontal  and  vertical  diffusion  displacemoits  are  written  as 
follows: 


L,  =  ^D^COSilnR) 


(8) 


=  )jD^SINQ.-KR) 


(9) 


L,  =  ^D,Ar(03  -  R) 


(10) 


where 
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=  X  direction  difliision  distance 
=  horizontal  diffusion  coefficient 
R  =  random  real  number  generator  between  0.0  and  1.0 
L^=Y  direction  diffusion  distance 
L,  =  Z  direction  diffusion  distance 
D,  -  vertical  diffusion  coefficient 

Subsequent  to  the  above  computations,  the  diffusion  distances  are  then  scaled 
and  added  to  tire  particle  update  equation  (Equation  1)  each  iteration. 

The  implementation  of  the  particle  tracking  model  into  the  CH3D  hydrody¬ 
namic  simulation  of  the  New  York  Bi^t  required  boundary  option  develop¬ 
ment  Boundary  options  were  developed  for  land  and  the  seaward  boundary. 
Two  options  are  available  for  land  boundaries  and  are  ^cified  by  means  of 
the  parameter  IBND.  When  IBND  is  set  to  zero,  a  particle  will  remain  on 
land  when  it  intersects  a  land  boundary.  In  this  case,  the  parameter  ISWC  is 
reset  from  zero,  an  active  particle,  to  one,  an  inactive  particle,  which  is 
removed  from  the  computation.  A  value  of  one  for  IBND  allows  a  particle  to 
slide  along  a  land  bouridary  and  remain  within  the  computational  grid.  Inflow 
and  outflow  qitions  are  available  for  seaward  or  water  boundaries.  A  particle 
ai^aching  a  seaward  boundary  passes  through  ^  boundary  and  is  removed 
from  the  computation  by  setting  tire  parameter  ISWC  to  one.  If  desired, 
inflow  boundaries  can  act  as  multiple  continuous  release  points.  This  feature 
is  controlled  by  the  parameters  ICON,  IFRC,  aixl  NPAR.  The  limitation  of 
this  option  is  that  it  only  works  for  point  continuous  releases,  which  means 
that  grid-wide  initial  particle  positicms  cannot  be  specified. 

The  particle  tracking  model  has  two  operational  options.  The  first  includes 
the  particle  tracking  model  as  subroutirres  within  CH3D.  In  this  mode,  uie 
particle  trajectory  computations  are  performed  in  concert  with  the  hydro- 
dynamic  simulation,  arid  the  display  of  particles  is  conducted  during  post¬ 
processing.  The  second  mode  performs  the  particle  tracking  computations 
within  a  stand-alone  graphics  program  that  allows  simultaneous  simulation  and 
animation  of  particle  movement  When  the  stand-alone  version  of  the  particle 
tracking  model  is  used,  CH3D  is  run  independently  and  the  time-invariant 
geometric  data  that  defines  the  grid  are  stored  along  with  temporally  averaged 
(e.g.,  1-hr  average)  velocity  components.  This  mode  has  the  advantage  that 
various  particle  release  points  can  be  easily  and  quickly  evaluated  using 
preprocessed  hydrodynamics.  Thus,  various  gaming  analyses  can  be  conducted 
such  as  deriving  impact  probabilities  utilizing  multiple  particle  release 
positions. 
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3  Stand-Alone  Particle 
Tracking  Model 


The  basic  idea  behind  the  stand-alone  panicle  tracking  model  and  grai^cal 
interface  is  that  CH3D  hydrodynamic  simulations  are  performed  independently 
where  time-invariant  geometric  informatirm  is  stored  along  with  temporally 
averaged  wind  speeds,  surface  drift  currents,  surface  elevatitHis,  arid  3D  veloc¬ 
ity  components  for  each  grid  cell  location  A  description  of  the  CH3D  hydro¬ 
dynamic  variables  can  be  found  in  Johnson  et  aL  (1^1). 

To  preserve  wind  speed  information  in  subroutme  CH3DWT,  two  additional 
arrays  have  been  added  to  CH3D  and  the  common  block  CH3D.INC.  The 
arrays  TXF  and  TYF  save  the  time  interpolated  values  of  the  dimensional 
Cartesian  wind  greeds.  The  arrays  TXFL  and  TYFL  contain  the  dimensional 
transformed  components  of  the  wind  speed.  The  wind  speed  data  saved  in  the 
TXFL  and  TYFL  arrays  are  used  by  the  floatable  force  balance  subroutine 
WFORCE  in  both  the  particle  tracking  subroutines  and  the  stand-alone  particle 
tracking  model  These  data  are  written  to  a  disk  file  during  the  CH3D  model 
run  by  subroutine  WQPM.  Stand-alone  particle  tracking  simulations  are  then 
performed  using  the  wirulows  interface  software  that  sets  up  and  executes  the 
particle  tracking  routine  PARS  and  the  graphical  output  routines. 

The  output  subroutine  WQPM  is  run  as  part  of  CH3D  during  a  hydrody¬ 
namic  simulation.  The  inputs  to  WQPM  (File  90)  are  the  averaging  interval 
NAVG  and  the  output  starting  iteration  FI^QS.  The  ou^t  from  WQPM 
(File  70)  are  (a)  time-invariant  transformation  data  that  geometrically  relate  the 
I^ysical  grid  system  to  the  transformed  computational  grid,  and  (b)  time- 
varying,  time-averaged  hydrodynamic  data  (File  71).  The  time-invariant  data 
are  ou4>ut  once  and  consist  of  rrondimensional  scaling  parameters  (UREF, 
XREF,  ZREF,  and  RB),  the  rrondimensional  hydrodynamic  simulation  time 
step  (DT),  spatial  transformation  derivatives  (XI,  X2,  Yl,  and  Y2),  metric 
tensors  (Gll,  G12,  G22,  and  GD),  and  grid  tow  and  column  reference  quanti¬ 
ties  (NROW,  IROW.  lUl,  IU2,  NCOL,  JCOL,  JVl,  and  JV2).  This  operation 
occurs  when  WQPM  is  first  called  by  CH3D,  or  when  the  time  iteration  coun¬ 
ter  IT  equals  FTWQS. 

During  each  CH3D  iteration  subsequent  to  ITWQS,  surface  velocities 
(USRF  and  VSRF)  are  ctnnputed  using  the  Bonan  surface  drift  subroutine 
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EK2,  and  averages  of  the  hydrodynamic  field  infoimaticm  (UA,  VA,  WA,  FIA, 
FJA,  TXFLA,  TYFLA.  and  AHSSA)  are  accumulated.  At  the  end  of  averag¬ 
ing  interval  NAVG.  the  particle  tracking  hydrodynamic  output  is  written  to 
File  71.  and  the  arrays  are  reset  to  zero. 

Extensive  testing  of  the  various  particle  tracking  options  has  shown  that  the 
small  amount  of  additioiud  computational  effort  requi^  by  the  higher  order 
interpolation  and  integration  is  warranted  when  compared  with  the  additional 
accuracy  gained.  As  a  consequence,  dre  particle  tracking  qrtions  selected  for 
use  in  die  stand-alone  particle  tracking  model  (PARS)  are  (a)  second-order 
trapezoidal  rule  time  integration,  (b)  biquadratic  spatial  interpolation  over  each 
horizontal  grid  ceU  widi  a  two  point  vertical  interpolation,  (c)  nonlinear  hori¬ 
zontal  convective  acceleration,  (d)  Ekman  surface  drift,  (e)  direct  wind  forcing 
oh  floatables.  (f)  constant  speed  particle  setding,  (g)  isotropic  particle  diffu¬ 
sion.  and  (h)  single  and  multiple  position,  instantaneous  arid  continuous  parti¬ 
cle  releases. 

A  graphical  user  interface  providing  a  point  and  click  operating  environ¬ 
ment  has  been  created  to  aid  die  user  in  the  operation  of  die  stand-alone  parti¬ 
cle  tracking  model.  The  user  may  start  the  interface  by  typing  in  "pt”  from  the 
directory  containing  die  grtqphical  user  interface  software.  This  interface  pro¬ 
vides  a  main  ctmtrol  panel  firom  which  the  user  may  perform  ir^xit  and  model 
execution  tasks  (Figure  1).  Prior  to  executitm  of  die  model,  two  input  tasks 
are  performed:  (a)  specification  of  paitide  release  sites,  and  (b)  input  of 
model  control  parameters.  While  specifying  a  sequence  of  release  sites,  the 
main  control  panel  allows  die  user  to  eidier  create  a  file  containing  new  release 
sites  and  types  (ie.,  single,  multiple,  and  continuous  particle  releases)  or  select 
a  pre-existing  file.  Whenever  die  user  elects  to  create  a  new  release  site  file,  a 
window  will  igipear,  containing  an  image  of  the  New  York  Bight  grid,  that 
allows  the  user  to  point  and  dick  the  location  of  release  sites  (Rgure  2).  A 
listing  of  existing  release  site  files  can  be  generated  via  pushbutton  on  die 
main  control  panel.  This  allows  the  user  to  sdea  and  gnqihically  view  a 
particular  release  site.  The  user  is  also  provided  with  the  crgiatnlity  of  ddeting 
a  release  site  file  with  the  release  site  manager  pushbutton. 

Spedfication  of  model  input  control  parameters  such  as  particle  density, 
drag  coeflidents,  and  wind-forcing-enabled  is  accomidished  throu^  interaction 
with  the  main  control  panel.  The  create  pushbutton  dlows  die  user  to  create  a 
new  control  parameter  file.  A  dialog  wiridow  with  a  listing  of  control  para¬ 
meters,  default  settings,  and  hdp  pushbuttons  will  appear  on  the  screen  (Fig¬ 
ure  3)  when  the  user  sdects  the  create  jmshbution.  The  he4>  pushbutton 
provides  detailed  information  on  each  parameter.  All  control  parameters  are 
forced  to  have  a  spedfied  range  of  values.  If  the  user  attempts  to  set  a  para¬ 
meter  to  a  value  outside  its  dlowable  range,  it  will  be  reset  to  its  default  value 
and  a  wanting  dialog  window  will  t^ipear.  The  user  is  also  provided  with  the 
ability  to  view  and  edit  an  existing  control  parameter  file  by  clicking  on  the 
control  parameter  file  view  pushbudotL  The  user  is  provided  the  option  to 
delete  an  existmg  control  parameter  file  by  clicking  on  the  control  parameter 
manager  pushbutton. 
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The  final  task  consists  of  creating  and  viewing  the  particle  trajectory  output 
file.  To  create  a  trajectory  file,  the  user  will  execute  the  model  by  clicking  on 
die  trajectory  file  create  pushbutton.  A  dialog  window  prompting  the  user  for 
the  particle  tracking  trajectory  file  will  appear.  Once  the  user  selects  a 
filename,  anodier  window  will  pop  vp  containing  two  file  selection  dialogs 
(Hgure  4).  The  user  must  then  select  a  release  site  and  parameter  control  file 
to  continue.  Once  this  has  been  done,  the  model  is  executed.  During  model 
execution,  a  window  is  displayed  that  provides  a  model  run  status  field  and  an 
abort  pushbutton. 

Subsequent  to  model  execution,  the  user  may  view  the  trajectory  output  by 
clicking  on  the  trajectory  file  view  pushbutton  on  the  main  control  panel.  This 
provides  a  file  selection  dialog  window  promi^g  the  user  to  select  a 
trajectory  file.  Next,  a  window  will  appear  that  allows  the  user  to  view  an 
animation  of  the  particle  trajectories  (Figure  S). 

The  particle  traddng  model’s  gn^cal  user  interface  software  was  written 
using  a  combination  of  X-windows  and  the  Silicon  Graphics’  graphical  pro¬ 
gramming  library  GL.  Its  use  is  limited  to  cominiter  platforms  supporting  both 
software  packages  such  as  the  Silicon  Grai^cs  woricstation  used  during  soft¬ 
ware  development 
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4  Model  Verification 


The  initial  set  of  test  simulations  consisted  of  comparing  paitide  trajectory 
results  with  the  analytical  solution  of  Lagrangian  residual  currents  in  a  two- 
dimensional  Oongitudinal-veitical)  dead-end  tidal  channel  of  omstant  width, 
depth,  and  density  Ganniello  1977).  The  CH3D  model  setup  included  constam 
density,  no  wind  stress,  a  linear  bottom  shear,  a  constant  21.7-cmVsec  eddy 
viscosity,  and  a  simple  12.42-hr  period  tide.  The  rectangular  grid  was  20  cells 
long,  10  cells  deq),  and  1  cell  wide  with  grid  spacings  of  Ax  =  3.5  km,  Az  = 
1.0  m,  and  Ay  =  3.S  km.  The  purpose  of  die  test  was  to  determine  if  the  parti¬ 
cle  tracking  schemes  predicted  the  correct  residual  transport  within  the  channel. 
A  single  particle  was  initially  released  in  each  of  10  vertical  layers  and  tracked 
for  1  tidal  cycle  subsequem  to  a  10-tidal  cycle  spirt-up.  The  residual  velocity 
was  then  computed  by  dividing  each  particle  displacement  from  the  initial  lon¬ 
gitudinal  position  by  die  tidal  period.  The  results  of  the  first-order,  bilinear 
particle  track  test  computation  are  presented  in  Figure  6,  in  which  die  pre¬ 
dicted  Lagrangian  residual  velocity  is  compared  with  the  analytic  solutitm 
(lanniello  1977)  at  a  locatitm  7  km  frimi  tidal  boundary.  Examination  of 
the  figure  shows  dose  agreement  between  model  and  analytic  solution. 

Drogue  trajectory  data  collected  by  the  U.S.  Environmental  Protection 
Agency  (EPA)  for  ^  years  1989  thrw^  1991  were  examined  to  determine  if 
surface  trajectories  data  were  available  within  die  New  York  Bi^  hydrody¬ 
namic  model  grid.’  Of  die  more  than  50  drogue  trajectory  files  provided  by 
EPA,  only  1  contained  drogue  positions  diat  redded  widiin  the  grid  for  an 
adequate  enough  time  to  perform  a  model  comparison.  This  surface  drogue 
entered  the  omnputational  grid  from  die  east  on  22  April  1991  and  exited  to 
the  southwest  around  22  May  1991.  During  diis  time  period,  die  surface 
drogue  moved  from  40.2  N,  70.4  W  to  38.5  N,  73.5  W. 

Wind  speed  and  directitm  data  during  the  simulation  period  were  obtained 
for  Kennedy  International  Airport,  Ambrose  and  two  National  Oceanic 
and  Atmos^ieric  Administration  (NOAA)  data  buoys.  Kennedy  International 
(KI)  and  Ambrose  Li^t  (AL)  lie  north  and  west  of  die  drogue  trajectory  with 
coordinates  40J  N  and  73.8  W.  NOAA  data  buoy  Station  12  is  located  at 


'  Pexsonal  Onmnunication.  1992.  D.  Pabst,  EPA  106-Mile  Drifdqg  Buoy  Stu^  Project 
Maiuser,  Washington,  DC. 
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38.8  N  and  74.6  W,  which  is  west  and  south  of  the  drogue  path.  NOAA  data 
buoy  Station  4  is  located  south  and  west  of  the  drogue  trajectory  at  coordinates 
39.0  N  and  70.0  W. 

To  determine  which  meteorological  data  set  or  comlnnations  of  data  sets 
stmuld  be  used  to  drive  the  particle  tracking  model,  hourly  data  for  paired 
locations  were  plotted  against  one  another  (Figures  7-12).  Examination  of 
these  figures  letuls  to  a  number  of  observations.  The  most  obvious  and 
troubling  observation  is  that  there  is  very  little  correlation  or  consistency 
between  stations.  For  example,  a  comparison  of  wind  speed  and  direction 
between  Ambrose  Light  and  Kennedy  International  (Figures  7  and  8)  shows 
that  the  wind  speed  at  Ambrose  Light  is  nonnally  greater  than  that  at 
Kennedy.  This  might  be  attributed  to  an  over-land  versus  over-water  wind 
effect,  and  there  is  considerable  scatter  in  the  observed  directions.  A  compari¬ 
son  of  an  off^re  and  nearshore  station,  NOAA  data  buoy  4  and  Ambrose 
light,  shows  that  the  nearshore  wind  speeds  are  usually  greater  than  that 
observed  offshore  (Figure  9).  Althou^  the  wind  directions  in  an  average 
sense  appear  to  be  similar,  there  is  considerable  scatter  in  the  hourly  measure¬ 
ments  figure  10).  Finally,  a  comparison  of  the  NOAA  data  buoys  4  and  12 
(Figures  11  and  12)  shows  that  the  offshore  wind  speeds  are  comparaUe,  but 
again  there  is  scatter  in  the  measured  directions. 

Wind  data  for  the  time  period  22  April  to  22  May  1991  were  available  at 
Kennedy  International,  Ambrose  Li^  and  NOAA  data  buoy  12.  Model 
predictions  utilizing  wind  speeds  and  directions  from  these  stations  did  not 
compare  well  with  the  observed  trajectory.  For  example,  die  predicted  particle 
trajectory  (dadied  line)  using  NOAA  data  Inioy  12  is  compared  widi  the 
drogue  observations  (solid  line)  in  Figure  13.  In  this  figure,  the  open  box 
denotes  the  starting  position  and  the  closed  circles  denote  the  end  position.  It 
is  seen  in  diis  figure  that  bodi  die  predicted  mean  direction  and  total  distance 
traveled  are  not  predicted  well.  The  inability  of  the  model  to  accurately  pre¬ 
dict  the  drogue  trajectory  using  any  of  the  available  individual  wind  data  can 
be  attributed  to  (a)  die  remoteness  of  the  wind  data  sets  stations  from  the  tra¬ 
jectory  positions,  and  (b)  the  lack  of  uniformity  of  the  grid  wide  wind  field  as 
exhibited  by  the  lack  of  correlation  between  wind  stations. 

In  an  attempt  to  generate  a  wind  field  that  is  more  representative  of  the 
large  scale  or  geostroidiic  wind  variability  during  the  simulation  period,  mean 
wind  speeds  and  directions  were  determined  from  the  offslmre  stations  after  all 
short  duration  or  localized  storm  events  were  removed.  In  other  words,  (xily 
mean  wind  and  directions  that  persisted  for  a  day  or  more  were  included  in  the 
constructed  wind  field.  A  repeat  of  the  particle  tracking  simulation  using  the 
constructed  wind  fidd  result^  in  a  predicted  particle  trajectory  that  more 
reasonably  mimicked  the  observed  data.  Table  1  lists  the  observed  and  pre¬ 
dicted  latitude  (N)  and  longitude  (W)  drogue  trajectory  positions,  respectively. 
In  this  table,  simulation  day  zero  corresponds  to  midnight  on  23  April  1991. 
Figures  14  and  15  provide  correlation  plots  of  the  data  presented  in  Table  1. 
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5  Conclusions  and 
Recommendations 


At  the  present  time,  the  implementation  of  the  panicle  tracking  model  is 
complete.  Additional  improvements  that  could  enhance  die  predictive  ct^iabil- 
ity  of  the  particle  tradting  model  are  (a)  a  "pufT  or  near  field  outfall  source 
algorithm,  (b)  statistical  reiHesentation  of  particle  size  or  density  for  differen¬ 
tial  settling  analyses,  (c)  anisotropic  panicle  diffusion  that  varies  with  the  qia- 
tial  scale  of  the  release,  and  (d)  near-bottom  bedload  tran^n  because  of 
currents  and  waves. 

With  respect  to  floatables  prediction,  given  that  the  offshore  transpon  of 
floatables  is  controlled  almost  entirely  by  the  local  winds,  it  is  imperative  that 
an  accurate  characterization  of  both  the  mean  wind  speed  and  direction  is  used 
to  drive  the  particle  tracking  model.  Consequently,  if  winds  remote  to  the  area 
of  interest  must  be  used  in  the  particle  tracking  simulation,  care  should  be 
taken  to  remove  localized  storms  or  varialnlity  from  the  wind  record.  If  possi- 
Ue,  barometric  pressure  (diarts  (ie.,  kinematic  analysis;  Cardone  1969; 

CardotK,  Greenwood,  and  Greenwood  1992)  should  be  used  in  conjunction 
with  the  point  wind  measuiemoits  to  ensure  that  the  long-term  variations  in 
mean  wind  speed  and  direction  ate  represented  in  tiie  q>plied  wind  field. 
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Figure  2.  Release  stte  specification  screen 


Figure  6.  Comparison  of  predicted  and  analytical  Lagrangian  velocities 
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Figure  8.  WirxJ  direction  comparison  between  Ambrose  Light  and  Kennedy  International 
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Figure  9.  Wirxl  speed  comparison  between  NOAA  Station  4  and  Ambrose  Light 


L  AL  (  Me  L ) 


400 .  00  -1 


300. 00 H 


o 

C 

o 

^  200 . 00  - 
u 

L 

Q  : 


■a 


100.00- 


*  *  * 


0.  00  ~|-Tn~n’n~r-r--r'rTy'i--’rTn-'"i  i  i  i  i"|  *rT’T7*rT*n-f  r*i"t  i  ~m -n-i  | 

0.00  100.00  200.00  300.00  400.  U.. 

Wind  Direction  at  Station  4  (  I'lo  t ) 


Figure  10.  Wind  direction  comparison  between  NOAA  Station  4  and  Ambrose  Light  (degrees 
from  north) 
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Figure  1 1 .  Wind  speed  comparison  between  NOAA  Stations  4  and  12 
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